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ABSTRACT 

The observed interstellar scintillation pattern of an intra-day variable radio 
source is influenced by its source structure. If the velocity of the interstellar 
medium responsible for the scattering is comparable to the earth's, the vector 
sum of these allows an observer to probe the scintillation pattern of a source in 
two dimensions and, in turn, to probe two-dimensional source structure on scales 
comparable to the angular scale of the scintillation pattern, typically ~ 10 //as for 
weak scattering. We review the theory on the extraction of an "image" from the 
scintillation properties of a source, and show how earth's orbital motion changes 
a source's observed scintillation properties during the course of a year. The 
imaging process, which we call Earth Orbit Synthesis, requires measurements of 
the statistical properties of the scintillations at epochs spread throughout the 
course of a year. 

Subject headings: galaxies: active — scattering — galaxies: structure 



1. Introduction 

Flux density variability on intra-day timescales has recently been detected in a few tens 
of extragalactic radio sources (Heeschen 1984; Witzel et al. 1986, Kedziora-Chudczer et al. 
2001). There is now overwhelming evidence that the principal mechanism responsible for 
intra-day variability (IDV) at centimeter wavelength is interstellar scintillation (ISS), rather 
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than intrinsic source variability (Kedziora-Chudczer et al. 1997; Dennett-Thorpe & de Bruyn 
2001a,b; Jauncey et al. 2001; Rickett et al. 2001; Jauncey & Macquart 2001). 

ISS is observed in compact (< 100 //as) radio sources due to the scattering of radiation 
in the Galactic interstellar medium (ISM) as it propagates toward earth. Inhomogeneities 
in the ISM impart phase structure on the wavefront, which leads in turn to a pattern of 
flux density deviations on the observer's plane. Flux density variability is observed in a 
scattered source if an observer moves relative to the source's scintillation pattern, ie. due 
to relative motion between the ISM, the source and the earth. The presence of such radio 
IDV implies extremely small source sizes (~ 5 — 100 //as), and high brightness temperatures 
(Tg ~ 10 13 — 10 15 K, e.g. Kedziora-Chudczer et al. 1997). This is well above the inverse 
Compton limit, and it is not clear whether such emission is explainable in terms of Doppler 
boosted synchrotron emission (Begelman, Rees & Sikora 1994). The properties — including 
the radio structures — of IDV sources are therefore of considerable interest. 

The scintillation pattern of an IDV source contains information on the structure of the 
source, which can be derived from the statistical properties of the scintillation pattern given 
the statistical properties of the scattering medium. This has already been demonstrated for 
sources undergoing interplanetary scintillation; the experiment by Cornwell, Anantharama- 
iah & Narayan (1989) is an elegant illustration of this technique. 

Some IDV sources exhibit large annual modulations in their variability timescales (Dennett- 
Thorpe & de Bruyn 2001a; Rickett et al. 2001; Jauncey & Macquart 2001). This is explained 
in terms of the earth's orbital motion; as the earth's velocity changes during its orbit about 
the sun, so does an earth-based observer's velocity relative to the ISM. Both the speed and 
direction change greatly if the speed of the scintillation pattern relative to the heliocentre, 
due to the intrinsic motion of the ISM, is comparable to the earth's speed. This effect is 
expected to be prevalent for IDV sources because the ISM velocity is comparable to earth's, 
but is not generally observed in the scintillations of pulsars because of their high intrinsic 
speeds. Their speeds, typically > 150 km/s are much greater than that of the earth or the 
ISM, and usually dominate the scintillation speed. 

The large variation in scintillation speed and direction over the course of a year makes 
it possible to observe two-dimensional structure in the scintillation pattern. Hence ISS can 
be used to "image" a source on 5 — 100 /ias scales. In this paper we show how to use changes 
in the earth's velocity around the sun to probe the scintillation pattern in two dimensions, 
and how to invert this information to obtain source structure. We refer to this technique 
as Earth Orbit Synthesis in light of its similarity to the conventional radio interferometric 
imaging technique of Earth Rotation Synthesis. 



- 3- 



As for the radio imaging technique of Earth Rotation Synthesis, Earth Orbit Synthesis 
assumes that the source structure remains relatively stable over the course of the obser- 
vations. This appears to be a good assumption for several IDV sources that have been 
well-studied over a number of years. B0917+624 appears to have been scintillating over the 
decade and a half that this source has been followed in detail (e.g. Jauncey & Macquart 
2001), while J1819+3845 has shown a clear annual cycle that repeats closely over at least 
the last two years (Dennett-Thorpe & de Bruyn 2001a), and PKS 1519-273 has shown IDV 
each time it has been observed in sufficient detail over the past decade (Kedziora-Chudczer 
et al. 2001). However, it should be remembered that this technique is not applicable to those 
sources that exhibit episodic IDV such as PKS 0405—385 (Kedziora-Chudczer et al. 2000). 

In §2 the apparent scintillation velocity of an extragalactic source is derived in terms 
of the velocities of the ISM and the earth. In §3 we relate the statistics of intensity and 
polarization fluctuations to source structure. A model for the scintillation velocity is com- 
bined with statistics of the intensity fluctuations to derive source structure in §4, and the 
technique is demonstrated with some specific examples. The conclusions are presented in §5. 



2. Earth's orbital motion 

The velocity of an earth-based observer across the scintillation pattern of a scattered 
source varies annually due to the earth's orbit about the sun. Here we derive the observed 
scintillation velocity in terms of the position of a distant scintillating source and the velocity 
of the interstellar medium relative to the heliocentre (see Fig. 1). 

We assume that the internal motions of the scattering material are negligible compared 
to the bulk motion of the material across the line of sight. We also assume that the scattering 
occurs predominantly in one layer of the ISM only, at a distance L from earth, so that motion 
of the scattering material can be characterized by a single velocity. This assumption appears 
to be valid for a few IDV radio sources whose scattering screens appear to be nearby and 
localized (e.g. Dennett-Thorpe & de Bruyn 2000), and is a good approximation for other 
IDV sources. 

The observed velocity of the line of sight to an extragalactic source through a scattering 
screen is given by 

Vobs = Vism - V0, (2-1) 

where vism = (fiSMxj ^ism^, wism 2 ) and are, respectively, the velocities of the ISM and 
the earth relative to the heliocentre. For extragalactic sources, the proper motions of the 
source across the line of sight can be ignored. The co-ordinate system is aligned so that e 2 
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points normal to the ecliptic plane and e x points along ecliptic longitude 0° in the ecliptic 
plane. 

We ignore the 0.01° inclination of the earth's orbit to the ecliptic, and approximate the 
earth's orbital motion as occurring in the ecliptic plane. The eccentricity of the earth's orbit 
(e = 0.0167) results in a negligible variation in orbital speed. We therefore approximate the 
earth's orbital velocity relative to the heliocentre as: 

= f0(sin27rt, — cos27rt, 0), (2-2) 

where vq = 29.78 km/s is the mean orbital speed of the earth, t is the time in years measured 
from the vernal equinox. 

As shown in Figure 1, and choosing the earth to be at the origin of the co-ordinate 
system we write location of the scintillation screen in terms of ecliptic co-ordinates a = 
a (cos <j) cos 9, sin cos 9, sin 9), where and 9 are, respectively, the ecliptic longitude and 
latitude of the scintillating source. For typical screen distances > 20pc, changes in a = a/a 
caused by orbital changes in the earth's position relative to the screen are negligible. 

The observed motion of the scintillation pattern across the line of sight is 

v scint = \Y\\ + v cos(27rf — <f>)]e\\ + [V± + v sin 9 sin(27rf — <f>)]e±, (2-3) 

where en = (— sin0, cos0, 0) and e_|_ = (— cos sin 9, — sin c/>sin 9, cos 9) are, respectively, 
the unit vectors parallel and perpendicular to the ecliptic along the line of sight to the 
scintillating source. The intrinsic speed of the ISM relative to the heliocentre and tangential 
to the line of sight is V\\ along ey and V± along e_|_, where 

V\\ = wisM y cos0 - vi SMx sin(j), (2-4) 
V± = vjsmzCosO -sin9(visM x cos(j) + visMySin(j)). (2-5) 

Equation (2-3) implies that the velocity of the scintillation pattern relative to the earth 
changes in both magnitude and direction on an annual cycle. 

3. Source structure 

The scintillation pattern of a scattered source is related to both its intrinsic structure and 
the inhomogeneities in the ISM. Source structure can be derived if the statistical properties 
of the scintillation pattern due to a point-source are known. In §3.1 we express for the 
scintillation pattern due to a point source in terms of the scattering medium's statistical 
properties. In §3.2 we show how the scintillation pattern is modified when the scintillating 
source has extended structure. 
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3.1. The scintillation pattern of a point source 

The scintillation pattern of a point source depends upon the phase fluctuations on the 
scattering screen which are characterized by the phase structure function 

D,(r) = ([0(r' + r)-0(rO] 2 ). (3-1) 

The angular brackets denote an ensemble average over all possible realisations of the scat- 
tering medium, and r is a two-dimensional vector on the scattering screen. 

For most lines of sight through the ISM the phase structure function reflects an under- 
lying power law spectrum of phase inhomogeneities between some inner and outer scales Iq 
and L Q respectively, and one has 

D ^)=(- L -Y 2 > (3-2) 

where (3 is close to 11/3 (Armstrong, Rickett & Spangler 1995), the value expected for 
Kolmogorov turbulence. The strength of the phase fluctuations is measured by rdiff, the 
length scale over which the root-mean-square phase difference on the phase screen is one 
radian. 

The strength of the scattering is determined by the ratio r F /r diff ; the quantity r F = 
\JcLj (2z/7r) is the Fresnel scale. The scattering is weak when phase irregularities on the 
scattering screen are small, when r^ff > r F and strong when they are large, when r^g < r F . 
The scattering strength varies as a function of frequency since both r diff oc z/ 2 /^" 2 ) and r F 
vary with frequency. It is convenient to define v t as the frequency at which the scattering 
changes from strong to weak. The scattering is weak at frequencies above is t , while strong 
scattering occurs at frequencies below u t . 

The qualitative properties of the scintillation depend heavily on the scattering strength, 
and hence on v t . In general, the value of v t toward a scintillating source is unknown, but can 
be determined from observations of its variability over a range of frequencies (e.g. Kedziora- 
Chudczer et al. 1997; Quirrenbach et al. 2000; Kedziora-Chudczer et al. 2001), or can be 
inferred from models of the ISM (Walker 1998). Along lines of sight out of the Galactic 
plane the transition frequency is in the range v t ~ 3 — 8 GHz (Walker 1998). 

A useful statistic of the intensity fluctuations of a point source due to scattering by the 
ISM is the intensity autocorrelation function 

C npt (r) = (5I pt (r' + r )5I pt (r')), (3-3) 

where 5I pt (r) = [I(r) — (/)]/(/) is the fractional deviation in the intensity at the point 
r on the observer's plane. This function describes the correlation between the intensity 
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fluctuations between points r' and r'+r on the observer's plane. An alternative but equivalent 
representation of the intensity fluctuations is provided by the power spectrum, Wn pt , which 
is related to the intensity autocorrelation function by a Fourier transform 

W IIpt = J d 2 rexp[tq • r]C /Jpt (r). (3-4) 

We adopt the power spectrum representation in the discussion below. 

The functional dependence of Wu pt (r) depends on the scattering strength. In the regime 
of weak scattering, one has (Salpeter 1967, Jokipii & Hollweg 1970, Narayan 1992) 

(3-2 /-, 

2 1 2 2 



W IIp M) = sin^-r^j, (3-5) 

with C = 11.2 for (3 = 11/3. The power spectrum peaks at q ~ l/r F , showing that the typical 
length scale of scintillations in the regime of weak scattering is r rp. An object undergoing 
weak scintillations would exhibit flux density variations on a timescale tp = rp/v sc i nt - 

In the regime of strong scattering the power spectrum of intensity fluctuations for (3 < 4 
is (Gochelashvily & Shishov 1975) 



W IIpt (q) 



-l 

ref 



2 ,- 3 ^-2)(4^)r W 2) q -2 ^_y~ ( , rF )6-/ 3 exp[-(,r ref )^ 2 ], q<r 

J^2 ^ (^2) r diS' r rcf ~ ? 

2^(/3-2)r(/3/2) n -2( n -\2-0 _-l < „ 



where r re f = rp/raiff is known as the refractive scale. Intensity fluctuations can occur predom- 
inantly on two timescales in this regime, corresponding to the two peaks in Wu pt (q). The 
outermost peak, corresponding to the effect of diffractive scintillation, occurs at q ~ l/r^is 
and is associated with fast-timescale, narrowband, large amplitude (up to 100% modulated) 
variability. 

The other peak in Wu pt (q) occurs at q ~ l/r re f, and is associated with refractive 
scintillation. Refractive scintillations occur on a longer timescale than weak scintillations, 
are broadband and have a typically low (~ 10%) modulation index. The timescale associated 
with refractive scintillation is t re f = r re f/i> S cint- This is a factor r-p/r^ > 1 slower than for 
weak scattering. 

Intensity fluctuations of IDV sources scintillating in the strong scattering regime are 
seen to occur on timescales ~ 1 — 10 days with amplitudes up to ~ 10% of the total source 
intensity (e.g. Kedziora-Chudczer et al. 2001). This is attributed to refractive scintillation 
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which is exhibited by any source with angular diameter 6 s < r Te f/L. No IDV source has yet 
been observed to exhibit diffractive scintillation. The stringent limitations on source size 
required for a source to exhibit diffractive scintillation, 6s < r^/L are harder to satisfy for 
IDV sources. For example, one requires #s < 6 /xas at v — 2 GHz to exhibit fully modulated 
diffractive scintillation scattering, where we assume v t = 5 GHz and D = 100 pc. These 
limitations are less stringent if D is smaller, as has been suggested for at least one source 
(Dennett-Thorpe & de Bruyn 2000). 

Expressions (3-5) and (3-6) are only valid when the statistical properties of the phase 
fluctuations on the scattering screen are isotropic, in the sense that the length scale rdiff 
does not vary with orientation on the screen. Anisotropy introduces a preferred direction 
in the scattering medium so that Cjj pt is no longer a function of r only, but depends on 
the orientation. Typical elongations of the scattering pattern in the ISM can exceed 2:1 
(see Chandran & Backer submitted for a full discussion). Anisotropic screen structure can 
therefore cause elongations in the scintillation pattern similar to apparent elongations caused 
by intrinsic source structure. One way to distinguish between source structure and screen 
anisotropy involves imaging the scattered source at low frequency, where the scattering is 
strong (ie. v < v t ) and the condition &s < r TC {/D is more likely to be satisfied. Then the 
shape of the scattered image is entirely dominated by the scattering and any anisotropy is 
readily apparent as an elongation of the scattering disk along the axis corresponding to small 
rdiff- This permits a direct measurement of the effects of screen anisotropy. Of course, the 
measurement can only be attempted at frequencies low enough that the scattering disk is 
sufficiently large to be resolved; this may be an important limitation in practice. 

We see in §4.1 that the effect of screen anisotropy is not important when the angular size 
of a source exceeds the angular scale of the scintillation pattern for a point source (ie. r F /D 
for weak scattering), and is readily distinguished from source structure by its characteristic 
structure in Wu pt . Thus, to first order, effects of screen anisotropy are diminished in many 
IDV sources. Hereafter we restrict our discussion to isotropic phase screens for the purposes 
of simplicity. 

It possible to determine which of equations (3-5) and (3-6) is applicable to a set of 
observations provided that the characteristics of the intensity fluctuations are known at 
several frequencies. For weak scattering one has tp oc z/ 1 / 2 , whereas t rc f oc v 2 " 2 (assuming 
(3 = 11/3) for strong scattering. The large difference in the timescales of the intensity 
variations allows a determination of whether the variability is occurring in the regime of 
weak or strong scattering. 

Many sources are observed to scintillate at frequencies that straddle the regime between 
weak and strong scattering. It is difficult to provide analytic expressions for Wu pt in this 
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transition regime, and instead, equation (3-5) may be used on the understanding that it is 



We now discuss the effect of source structure on the observed scintillation pattern of an 
extended source. Direct measurements of the scattered wavefield from an extended source 
may be inverted to completely determine the source structure. The intensity scintillation 
pattern of a source also contains information on source structure. However, analysis of the 
intensity fluctuations alone only provides partial information on source brightness distribu- 
tion. 



The wavefield of a scattered source is directly related to its structure. The mutual 
coherence function is the simplest nontrivial moment of the wavefield. If the wavefield 
measured by a receiver at the point r on the observer's plane is u(z;r), the normalized 
two-point correlation function is then 



where a, a' = [x,y] corresponds to the component of the wavefield oscillating along the x 
and y axes orthogonal to the direction of propagation, along the 2-axis, and J is the mean 
source total intensity. The normalized Stokes visibilities are defined by 



T I (z;r) = (r xa .(z;r) +r ro (z;r))/2 F Q (z; r) = {Y xx (z\ r) - r yy (z; r))/2 
rV(z;r) = (Y xy (z-v) + Y yx (z-v))/2 T v (z; r) = i(T xy (z; r) — T yx (z; r))/2. (3-8) 



It is assumed that the source is spatially incoherent and that the wavefront from the 
IDV source incident upon the scattering screen is planar, which is a very good approximation 
for extragalactic sources. (The planar approximation is not valid for Galactic objects, but it 
is possible to make a simple correction for the curvature of the wavefront in this case (e.g., 
Goodman & Narayan 1989).) The mutual coherence function measured by an observer is 



only approximately correct, or Wji pt may be derived numerically from scattering simulations. 



3.2. The scintillation pattern of an extended source 



3.2.1. The wavefield 




(3-7) 
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then related to the mutual coherence function of the source T x (0; r), X = [I, Q, U, V] by 



( lMz;r) \ 

r Q(^; r ) 



= e 



-A»(r)/2 



I r 7 (0;r) \ 

Tq(0; r) cos 2(f) V + Tu(0; r) sin 2<j) V 
-r Q (0; r) sin 2<j) V + Tu{0; r) cos 20 y 

V r ^(0;r) / 



(3-9) 



where <pv — A 2 RM is the mean phase change induced by Faraday rotation along the ray 
path, and we assume no inhomogeneity in the Faraday rotation transverse to the line of 
sight. 

The normalized mutual coherence of the source, Tx{0; r) is related to its angular bright- 
ness distribution normalized by the mean intensity via a Fourier transform 



r x (o ;r ) = / d 2 6b x (e)exp 



-2iri- 



6 ■ r 



A 



(3-10) 



The behaviour of the mutual coherence function depends on the angular size of the source 
relative to the angular scale of the scintillation pattern as seen by an observer on the earth. 
Equations (3-9) and (3-10) imply the behaviour of T(z; r) is dominated by the source angular 
structure if its angular size exceeds that of the scintillation pattern seen by an observer, 
9 S > k/r diS . In this instance the source is said to be resolved by the scattering medium. 

Conversely, the mutual coherence function of a source with angular size 6s < k/r^B is 
dominated by the scattering, and T(z; r) ~ T(z; 0) exp[— D ( f > (r)/2]. It is difficult to derive 
source structure when the source is unresolved by the scattering medium because the scat- 
tering properties dominate the mutual coherence function. This limitation is important in 
the regime of strong scattering where k/r&s is may be large relative to the size of an IDV 
source. 

The quantities Tx(r) can be measured over short coherent integrations (see, e.g. Good- 
man & Narayan 1989); typically two or three scintles are sufficient to obtain a good estimate. 



3.2.2. Intensity fluctuations 

Measurement of the intensity fluctuations contains partial information on the structure 
of a source. In this section we derive the relationship between the brightness distribution of 
a scattered source and its intensity fluctuations. The pattern of intensity fluctuations from 
an extended source is derived by dividing the source into small elements and computing 
the individual intensity scintillation pattern from each element (again assuming a spatially 
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incoherent source). The scintillation pattern of the extended source is therefore the sum of 
the scintillation patterns from all the elements of the source. 

Figure 2 shows that the wavefront incident on the screen at an angle 9 to the line of sight, 
and heading towards the point r + LO, intersects a thin scattering screen at the same point 
as a normally incident wavefront, heading towards r. Each of these wavefronts experiences 
identical scattering, and the amplification caused by the scattering screen is thus identical 
for these two wavefronts. It follows that scintillation pattern of an extended source is (Little 
& Hewish 1966), 



where SSx = (Sx — (Sx))/{Sx) is the fractional flux density deviation in Stokes parameter 
X = [I, Q, U, V], and b x (6) is the angular brightness distribution of the source as observed 
at the scattering screen, normalized such that J b I (6)d 2 6 = 1 and J b Y (6)d 2 6 = Y/I , Y = 
[Q, U, V\. The brightness distributions 6q and by at the screen differ from their distributions 
at the source if Galactic Faraday rotation is important. The distributions at the source are 
found by determining the Faraday rotation along the line of sight to the source. A magnetized 
ISM may also induce fluctuations in the amount of radiation transferred between Stokes Q, 
U and V, however the fluctuations in the degree of birefringence are of order fl e /uj p ~ 10~ 8 
smaller than the phase fluctuations considered here, and are probably unimportant for the 
scattering considered here (see Macquart & Melrose 2000). 

The flux density deviations observed in an extended source are connected to the scin- 
tillation statistics of a point source by 



The Fourier-transformed equivalent of equation (3-12) illustrates clearly the connection be- 
tween source structure and the scintillation statistics via the power spectrum of intensity 
fluctuations, 




(3-11) 




(5S x (r')5S Y (r + r')) 

[ d 2 ed 2 e , b x (e)b Y (e')(5l pt (r' - Le)5l pt (r' + r - Le')). 



(3-12) 



W XY (q) = C XY (q) 



(S*(q)5£(q)> 
hBx^By^Wu^)], 



(3-13) 



where B x (0) = b x (6/L), the Fourier transform of a function is denoted with a tilde, and we 
have made use of equation (3-4). The quantity B x {0) is identical to the visibility T x (0; r) of 
the previous section. Equation (3-13) may be inverted to derive information on the structure 
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of a scintillating source if the scattering properties of the ISM, embodied in Wu pt (q), are 
known. 

Equation (3-13) may be interpreted as follows. The quantity W IIpt (q) is the power 
spectrum of intensity fluctuations in the scintillation pattern that would be observed for a 
point source. The quantity Sx(q)-By(q) is the power spectrum of intensity structure in the 
source. The resultant power spectrum of intensity fluctuations, Cxy(q), is that of a point 
source if the angular structure contained in source is smaller than the scale of the scintillation 
pattern for a point source (ie. the source is unresolved by the ISM). On the other hand, the 
observed scintillation pattern is dominated by source angular structure if the source is larger 
than the scale of the scintillation pattern for a point source. For example, scintillations in 
the weak scattering regime contain power on length scales r > r F (equivalent to q < l/r F ). 
Source structure has no measurable effect if the source angular size is less than the angular 
scale of the scintillation pattern (ie. when 9 S < r F /D). 

If the angular size of the source is smaller than the scale size of the scintillation pattern 
for a point source, the term in square brackets in (3-13) is approximated by Wn pt (q). This 
is because M / 7/ pt (q) decreases much faster with q than _B/(q)i?/ (q). On the other hand, if 
the source angular size exceeds that of the scintillation pattern for a point source, we can 
use the approximation Wn p t ~ 2.8 g 4 ~^r F /r^ 2 for small q to obtain (see equation (3-12)) 



Sources sufficiently large to be resolved by the scattering medium exhibit smaller and 
slower intensity fluctuations relative to unresolved sources. The intensity modulation in- 
dex, Cxy(O), of an resolved source is therefore smaller than that of an unresolved source. 
This is evident from equation (3-13). The amplitude of Cxy(O) decreases if the terms due 
to the source cut off Cxy(q) before Wjj pt reaches its maximum value (e.g. at q ~ l/r F for 
weak scattering). 



The earth's orbital motion allows the scintillation pattern of a scattered source to be 
mapped in two dimensions. This is because the velocity of an earth-based observer relative 
to the ISM changes as the earth's orbital velocity changes during the course of a year. An 
observer moving at a velocity v relative to the ISM samples the scintillation pattern along the 
direction parallel to v. It is possible to measure the correlations T x (r = vt) or Cxy(r = vt) 
for observations a time t apart. The quantity T x (r) can be measured over short coherent 




(3-14) 



4. Earth Orbit Synthesis 
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integrations (see, e.g. Goodman & Narayan 1989). A good approximation to the quantity 
Cxy( r ) can typically be derived from observations of > 5 scintles, typically requiring one to 
two days flux density monitoring for typical IDV sources. 

Over the course of the year observations at different epochs sample the scintillation 
pattern in different directions, as the velocity of the earth changes around the solar circle. In 
particular, if the velocity of the ISM is comparable to the earth's, then a strong "annual cycle" 
effect will appear in the observed IDV pattern, as has been observed already in J1819+3845 
(Dennett-Thorpe & de Bruyn 2001a) and 0917+624 (Jauncey & Macquart 2001; Rickett et 
al. 2001). This enables measurement of Tx(r) or Cxy(r) for a large range of orientations of 
r, and the information may be inverted to obtain source structure provided that it is stable 
on a timescale long enough for it to be mapped. The main limitation to the technique is that 
one does not always know v because the velocity of the ISM is not generally known. Thus 
one must solve simultaneously for v and the source structure, or measure v directly from 
observations of the time delay between the scintillation pattern as observed at two widely 
separated (> 3000 km) receivers as the scintillation pattern passes across the earth (Jauncey 
et al. 2001, Dennett-Thorpe & de Bruyn 2002). 

We demonstrate the technique with several examples, concentrating on Earth Orbit 
Synthesis using intensity fluctuations. 



4.1. An example using Intensity Fluctuations 

An intensity image of the structure in a scintillating source requires measurement of 
the autocorrelation function C//(r). It is not possible to measure the spatial structure of 
the scintillation pattern of a source, C//(r), directly. The scale of the scintillation pattern 
is usually larger than that of earth so that it is impossible to measure Cn(r) on sufficiently 
long baselines r. This problem is overcome by measuring the correlations of the scintillation 
signal as a function of time, rather than baseline. One measures the temporal correlation 
function 

(/(W + 0> 

The correlation between two scintillation signals measured a time t apart at the same receiver 
may be derived from the correlation in the scintillation pattern measured by two receivers 
a distance r apart. The relation between the temporal and spatial intensity autocorrelation 
functions is thus 

C„,(t) = C„(r = vt) = p-Jj^ J d\e-">' [S(q)S;*(q)W> /pl (q)] . (4-2) 
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An important point in connection with (4-2) is that the modulation depth of the scintillation 
pattern, Cj/(0) is independent of the velocity of the scattering medium relative to earth. For 
instance, the depth of the intensity modulations does not depend on whether the scintillation 
pattern moves parallel or orthogonal to the long axis of an elongated source. However, we 
see that the characteristic timescale of the scintillations depends strongly on the direction 
of motion if v varies during the course of a year. 

The effect of source structure is particularly simple when a source undergoes weak 
scattering and its angular size is larger than that of the scintillation pattern of a point 
source (ie. $s > r-p/D). Equation (3-14) is applicable in this case, and we assume /3 — 4 
instead of the usual Kolmogorov value of (5 = 11/3 to simplify the algebra and derive analytic 
results. We obtain 

Cn(r) = ^| / d\e^ ^(q)! 2 , L6 S > r F . (4-3) 
r difr- L J 

This expression, and the assumptions made in deriving it, are employed in the following 
examples. 



4-1.1. Elongated source 



We illustrate Earth Orbit Synthesis by examining the combined effect of source structure 
and velocity variations on a source with elliptical brightness distribution: 



He) = 



exp 



Ox 

a. 



a, 



(4-4) 



where a x and a y are chosen to be the source angular scales parallel to the en and e_|_ axes 
respectively. This choice is made for convenience only, and it is simple to re-orient the source 
at an angle £ to the axes by making the transformation (x, y) — > (x cos £ + y sin £, —x sin £ + 
ycos£). Combining equations (4-3) and (4-4) yields the following intensity autocorrelation 
function 



Cn{v) = 



2.84 



2Tra x a y L 2 rl m 



exp 
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1/ 
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La X:V > r F 



(4-5) 



It is possible to determine the temporal decorrelation of the scintillation signal, for a given 
scintillation speed and direction. In the following we use V\\ = 15 km/s, V± = 10 km/s and 
we assume the source has ecliptic co-ordinates 6 = 30° and = 0. The variation in the 
scintillation speed over the course of a year is shown in Fig. 3. Figures 4(a) and 5(a) show 
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how the intensity temporal autocorrelation function Cu(t) varies throughout the year for an 
elliptical source with these screen parameters. 

The effect of source structure is evident in the timescale of variability. It is convenient 
to define the decorrelation timescale, tdc as the time in which the intensity correlation falls 
to 1/e of its maximum value. The decorrelation timescale for the elliptical source described 
by equation (4-5) is 

-1/2 

(4-6) 



£dc 



(v SC i n t • ey) 2 (v scint • e ± f 



2L 2 a 2 x 2L 2 a 2 y 



Changes in the magnitude of v scint cause changes in the scintillation timescale even if the 
source is circularly symmetric, whereas changes in the direction of v scint only result in changes 
in tiss if the source is elongated, when a x ^ a y . 

Figures 4(b) and 5(b) plot the variation in decorrelation timescale as a function of 
epoch for a specific scintillation velocity for the elliptical source with two different source 
parameters, La x = 10 9 m and La y = 10 9 m and 5 x 10 9 m respectively. Fig. 4(b) shows the 
expected scintillation timescale for a circularly symmetric source, and Fig. 5(b) the timescale 
for an elongated source. The peak in the scintillation timescale at T « 0.9 yr is common to 
both plots and occurs because the scintillation speed reaches a minimum. The emergence 
of another peak in Fig. 5(b) occurs due to the presence of structure in the source, here an 
elongation along the y-axis. The peak occurs at the time of the year that the scintillation 
direction is tangential to the long axis of the source. 



4-1-2. Two-component source 



A source with two components shows richer structure in the statistics of the intensity 
variations. We consider two gaussian patches of emission, the central peak having an intensity 
cto/o (0 < a < 1) and the other, located at an angular separation of Oq, having an intensity 
(1 — ao)I . The two subsources have identical angular scales a and we again assume that the 
angular diameters of the subsources exceed the angular diameter of the scintillation pattern 
due to a point source, so that equation (3-14) applies in the regime of weak scattering. The 
source angular brightness distribution is 
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(4-7) 



for which the resulting spatial intensity autocorrelation function is 



Cn(r) 
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(4-8) 
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Figure 6 plots the spatial intensity autocorrelation function of this source. This function 
has three peaks. The lightcurve of a two-component source therefore exhibits interesting 
behaviour when the scintillation direction is parallel to the line joining the two substructures. 
An observer then measures the spatial autocorrelation function along a line that joins the 
three peaks in the autocorrelation function. 

A specific example of such behaviour is illustrated in Figure 7, which demonstrates how 
the temporal autocorrelation function of a two component source changes throughout the 
course of a year. At T pa 0.0 — 0.2 yr the scintillation velocity is parallel to the direction of the 
two separated peaks, and a second bump is visible in the temporal intensity autocorrelation 
function. 



4-1.3. Polarization reference mapping 

Only limited information about source structure is available from Earth Revolution 
Synthesis using intensity fluctuations. It is impossible to determine bi(0) directly from 
C/j(r). Only the correlation of the angular brightness distribution, (6/(0 + 0') 6/(0')), can be 
derived. 

However, Earth Orbit Synthesis is useful even when the brightness distribution bj(8) 
cannot be ascertained exactly. IDV sources often exhibit fluctuations in Stokes Q, U and 
V that are dissimilar to those observed in /. This occurs if the polarized structure of the 
source differs from its structure in total intensity. This substructure can be mapped relative 
to the total intensity structure, using the correlations Cjq(t), C w (y) and C IV (r) to find the 
quantities (b I (0' + 0)5q(0')), etc.. It is thus possible to map the brightness distribution of a 
source in one Stokes parameter relative to another Stokes parameter, as has been shown for 
both the linear and circular polarized components in PKS 1519—273 (Macquart et al. 2000). 

5. Conclusion 

Measurement of the scintillation characteristics of IDV sources allows their structure 
to be investigated with microarcsecond resolution. The radio structures of IDV sources are 
currently of interest because of their high brightness temperatures. Earth Orbit Synthesis 
has the potential to determine the origin of these compact structures; it can determine, for 
example, whether the observed emission emanates from the central cores of these objects or 
from the bases of jets. 

Short timescale polarization variability in IDV sources indicates the presence of polar- 
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ization structure on similarly fine scales and — in many cases — on even smaller scales 
than the structure in total intensity. The advantage of Earth Orbit Synthesis over tradi- 
tional polarimetric imaging is its high resolution. Such a high resolution probe is less subject 
to beam-depolarization. Indeed, this mapping technique is already discovering remarkable 
detail about the polarized sub-structure of IDV sources (e.g. Macquart et al. 2000). 

The observed two-dimensional scintillation pattern of a scattered radio source is highly 
influenced by source structure. This pattern can be measured in two dimensions if the 
apparent velocity of the scintillation across the line of sight to a scattered source changes; 
changes in direction allow an observer to probe source structure along different directions. 
Direction changes occur when the velocity of the scattering material is comparable to the 
earth's velocity, and the vector sum of these two velocities changes during the course of a 
year. 

The derivation of source structure requires measurement of the scintillation fluctuations 
at several epochs over the course of a year. This is achieved by measuring the temporal 
decorrelation of the intensity {I(t)I(t + t')). Many IDV sources also exhibit polarization 
variability, which is often not well correlated with the total intensity. This occurs when 
the source structure in polarized and total intensity emission differs. Polarized substructure 
information is available by measuring the temporal cross-correlations (X(t)Y(t+t')) between 
Stokes parameters X and Y . It is possible to measure the polarized substructure of a source 
with reference to the intensity emission to high precision. 

An advantage of the Earth Orbit Synthesis technique is the capability to compare source 
structure between frequencies. For example, one could compare the angular separation 
between source components as a function to yield information relating to opacity changes 
with frequency. Analysis of the data on PKS 1519—273 (Macquart et al. 2000) demonstrates 
not only the power to compare source structure at two frequencies, but even to compare the 
relative source properties in different Stokes parameters. 

It should be emphasised that the use of Earth Orbit Synthesis as an imaging technique 
is subject to limitations. Firstly, the scintillation process only probes source structure over a 
limited angular size range. Consequently, while Earth Orbit Synthesis probes microarcsecond 
structures, it provides no information on the relation between microarcsecond and larger scale 
structure. Secondly, the use of autocorrelation functions in the imaging process is a necessary 
but severe limitation to any imaging that can be achieved. This technique yields the power 
spectrum of the source brightness distribution and, as such discards phase information that 
is present in normal rotational synthesis. The stochastic nature of the scintillation process 
forces us to use ensemble average quantities such as autocorrelation functions. 
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The angular resolution achievable using scintillation-based imaging depends on the 
properties of the scattering medium. In the regime of weak scattering, generally appli- 
cable at frequencies v > 5 GHz, the intensity scintillations are sensitive to angular scales 
6 ~ v-p/D oc v~ x l 2 . The scattering is strong at frequencies v < 5 GHz. The refractive 
scintillation exhibited by IDV sources in this regime is less sensitive to small angular scales 
because the angular resolution achievable is 6 ~ r TC f/D which is large, and scales as z/~ 2 - 2 
for a Kolomogorov spectrum of inhomogeneities. Strong scattering is typically sensitive to 
sizes > 0.1 mas, and hence can be found in many more flat-spectrum sources, particularly 
at low Galactic latitudes. Thus, while less useful competing against VLBI for milliarcsecond 
imaging, strong scattering here remains an excellent probe of the ISM close to the Galactic 
plane. 

A program to use Earth Orbit Synthesis to map structure in IDV sources is currently 
underway using the Australia Telescope Compact Array (Bignall et al. 2001). 

The authors thank Barney Rickett, Simon Johnston and Lucyna Kedziora-Chudczer 
for valuable discussions and comments on the manuscript. JPM is particularly indebted 
to Don Melrose, from whose Festschrift the idea for this paper originated. The Australia 
Telescope is funded by the Commonwealth Government for operation as a national facility 
by the CSIRO. 
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Fig. 1. — The scattering model. The earth-source line of sight moves relative to a scattering 
screen, fixed at a distance L from earth. 
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Fig. 2. — The geometry of the scattering assumed in calculating the effect of a source's 
structure on its scintillation pattern. 
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Fig. 3. — The scintillation speed of a source located at ecliptic co-ordinates (9, (f>) = (30°, 0°) 
with an intrinsic ISM velocity of (V\\, V±) = (15, 10) km/s. 
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300 T (years) 

Fig. 4. — The observed intensity autocorrelation function as a function of time and a epoch 
for an elliptical source whose angular size exceeds scattering disk and with a brightness 
distribution given by equation (4-4) with La x = 10 9 m and La y = 10 9 m. Panel (a) shows 
the evolution of the temporal autocorellation function, Cut throughout the year in the limit 
(3 — > 4, while (b) shows the evolution of the decorrelation timescale. 




300 T (years) 

Fig. 5. — The observed intensity autocorrelation function as a function of time and a epoch 
for a source whose angular size exceeds scattering disk and with a brightness distribution 
given by equation (4-4) with La x = 10 9 m and La y = 5 x 10 9 m. Panel (a) shows the 
evolution of the temporal autocorellation function, Cn t throughout the year (in the limit 
(3 — > 4), while (b) shows the evolution of the decorrelation timescale. 
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Fig. 6. — The spatial intensity autocorrelation function of a source whose angular size exceeds 
that of the scattering disk and with a brightness distribution given by equation (4-7) with 
La = 10 9 m and a source separation of L(9o x , 0y ) = (3, 0) x 10 9 m, a = 0.5 and assuming 
0-4. 
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Fig. 7. — The observed intensity autocorrelation function as a function of time and a epoch 
for a source whose angular size exceeds scattering disk and with a brightness distribution 
given by equation (4-4) with a x = 4 x 10 9 , a y = 10 9 and ao = 0.5. Panel (a) shows the 
evolution of the temporal autocorellation function, Cut throughout the year (/3 — ► 4), while 
(b) shows the evolution of the decorrelation timescale. 



